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10 years from now reactor neutrino experiments will attempt to determine which neutrino mass 
eigenstate is the most massive. In this letter we present the results of more than four million de- 
tailed simulations of such experiments, studying the dependence of the probability of successfully 
determining the mass hierarchy upon the analysis method, the neutrino mass matrix parameters, 
reactor flux models and, in particular, combinations of baselines. We show that a recently re- 
ported spurious dependence of the data analysis upon the high energy tail of the reactor spectrum 
can be removed by using a weighted Fourier transform. We determine the optimal baselines and 
corresponding detector locations. For most values of the CP-violating, leptonic Dirac phase S, a 
degeneracy prevents NO^A and T2K from determining either S or the hierarchy. We determine the 
confidence with which a reactor experiment can determine the hierarchy, breaking the degeneracy. 



Last year the Daya Bay [1, 2] and RENO [3] ex- 
periments demonstrated that the neutrino mass ma- 
trix mixing angle #13 is nonzero and as much as an 
order of magnitude larger than had been suspected 
several years ago. With the discovery that #13 is 
nonzero, at least three qualitative questions remain 
to be answered in the standard model plus three 
massive neutrinos. First, it is not known whether 
the second or third neutrino flavor is the most mas- 
sive. This choice is known as the neutrino mass hier- 
archy. Second, it is not known whether the leptonic 
sector has a non-zero CP-violating Dirac phase S. 
Third, the most precise determination of the mixing 
angle, O23, by the MINOS collaboration is 

sin 2 (2<9 23 ) = 0.96 ±0.04. (1) 

This equation has eight distinct solutions, known as 
the octants of #23 ■ These solutions can be divided 
into two families, corresponding to #23 modulo 90° 
greater than or less than 45°. 

The large value of #13 means that an appearance 
experiment, searching for v e 's in a beam of accelera- 
tor v^'s, may be sensitive to CP violation in the lep- 
tonic sector. The trouble is that, at the baselines of 
the T2K experiment in Japan and the NO^A exper- 
iment in the US, the CP phase 6 is degenerate with 
the mass hierarchy, the octant of 623 and the pre- 
cise value of O13. Following the analysis in Ref. [5], 
by combining data from the appearance of v e 's in a 
f M beam with the appearance of F e 's in a V ^ beam 
in its second three year run, given sufficient funding 
NO^A may be able to separate the octant and #13, 
which determine the total v e + v e appearance, from 



the hierarchy and S, which roughly determine the 
difference v e — V e . However for most values of 5 even 
the combination of T2K and NO^A will not be able 
to break the degeneracy between S and the hierarchy, 
and so will not be able to determine either. 

Fortunately the large value of #13 also implies that 
1-3 neutrino oscillations, those with amplitudes pro- 
portional to sin 2 (2#i3), are large enough to be ob- 
served in reactor neutrino experiments at medium 
baselines (40-60 km). These oscillations arc almost 
periodic in the inverse energy E, but their small 
aperiodicity may be used to determine the neutrino 
mass hierarchy [6] thus breaking the degeneracy at 
T2K and NO^A and so allowing a determination of 
the CP-violating Dirac phase. More precisely, com- 
bined with a hierarchy determination, the NO^A ap- 
pearance mode can only determine sin(<5) because 
the transformation 5 — > tt — S is degenerate with 
a slight shift in 623 or 813. However the effective 
mass difference determined by a reactor experiment, 
which will be discussed below, differs from the in- 
dependent atmospheric effective mass [7] determined 
by the disappearance channel of experiments such as 
MINOS and NOi/A. This difference may yield a la- 
determination of cos(i5). Such reactor experiments 
are not only possible but indeed will be performed 
within the next decade [8, 9]. 

In this letter we present the main results of a se- 
ries of simulations of medium baseline reactor ex- 
periments. We will determine the reliability of the 
determinations of the hierarchy in such experiments 
and will also find the optimal baselines for their de- 
tectors. In most of our simulations we use the value 
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of AM 2 2 ! from Rcf. [10], sin 2 (26»i 2 ) from Rcf. [11], 
|AM| 2 | determined by combining v and v mass dif- 
ferences from Ref. [4] and sin 2 (26*13) from [2] 

AM 2 \ = 7.59 x 1CT 5 cV 2 , sin 2 (2<9 12 ) = 0.857 
|AM 3 2 2 | = 2.41 x 10~ 3 eV 2 , sin 2 (2<9 13 ) = 0.089. 

Since we hold these parameters fixed, our neural net- 
work is trained with the same parameters which are 
used to test it, which artificially inflates the chance p 
of correctly determining the hierarchy. We have also 
systematically studied the effects of changing these 
parameters, as will be described below. 

We consider reactor experiments in which F e 's are 
detected via inverse f3 decay by a detector with a 20 
kton target mass consisting of 10% hydrogen. This 
is the mass of the proposed Daya Bay II detector, 
but it is four times the mass of the proposed de- 
tector for RENO 50, and so to interpret the results 
below for RENO 50 one needs to multiply all times 
by a factor of four. We assume a detector resolu- 
tion of 3%/ y/E(MeV). We ignore backgrounds and 
assume a perfect detector energy response. As a re- 
sult of these approximations our results are overly 
optimistic. 

We have performed 3 and 6 year simulations us- 
ing several reactor flux models with various cutoffs 
using fluxes arising from over 100 combinations of 
baselines. For each combination, we have have sim- 
ulated 5,000 experimental runs with each hierarchy. 
The data analysis methods described below are then 
used to attempt to determine the hierarchy. We re- 
port the percentage p of experiments for which the 
hierarchy is determined correctly. 

The determination of the neutrino mass hierarchy 
from 1-3 oscillations in the V e spectrum proceeds as 
follows. At energies E/McV greater than L/12 km, 
deviations from periodicity in l/E are too small to 
be measured and the wavenumber k determines [7] 

AM 2 ff = cos 2 (M| AM31 1 +sin 2 (6- 12 )|AM 2 2 |. (2) 

At low energies the deviation from periodicity is 
large and it determines various combinations of 
I AM. 2 ! I and |AAff 2 | given in Ref. [12], for example 
the energy of the 16th oscillation peak is propor- 
tional to IAM3J. To determine the mass hierarchy 
one needs to combine two distinct combinations of 
I AMgj I and | AM| 2 1 , thus one must combine the high 
and low energy parts of the spectrum. 



The most studied algorithm which determines the 
hierarchy given a reactor V e spectrum is that of 
Ref. [13]. One first finds the Fourier transform of the 
measured spectrum, as suggested in Ref. [14]. 1-3 os- 
cillations have a wavenumber k of about |AMf 2 | and 
so the peak structure of the transform at k <~ | AM| 2 1 
is sensitive to these oscillations and thus the hierar- 
chy. In Rcfs. [13, 15] it was shown that the heights 
of the peaks can be combined into two real num- 
bers RL and PV such that RL + PV is positive 
if and only if the hierarchy is normal. In Ref. [12] 
two more observables were added, one mixing infor- 
mation from the Fourier sine and cosine transforms 
and one using a nonlinear Fourier transform with the 
same aperiodicity as the 1-3 oscillations. A neural 
network finds the combination of these observables 
which best determines the hierarchy. 

An obstruction to this analysis has been described 
in Ref. [16], which observed that RL + PV is very 
sensitive to the choice of model of the reactor neu- 
trino flux and to variations of AM 2 ff smaller than 
the precision of its determination by MINOS [17]. 
To eliminate this spurious dependence, we use a 
weighted Fourier transform in which higher ener- 
gies are weighted less heavily, providing a soft cutoff. 
Fig. 1 shows the average values of RL + PV in sim- 
ulated data using a normal and a weighted Fourier 
transform. The normal Fourier transform leads to 
the fluctuating dependence found in Ref. [16], but 
these fluctuations are much smaller with a weighted 
transform. A steeper weight would further reduce 
the amplitude of the fluctuations, but would also re- 
duce the average value of RL + PV, weakening the 
hierarchy dependence. 

Thus we find that the spurious dependence ob- 
served in Ref. [16] can be eliminated via a sim- 
ple modification of the analysis. We will now dis- 
cuss a more serious challenge which afflicts otherwise 
promising sites for such an experiment. 

As a determination of the hierarchy requires a 
medium baseline, the neutrino flux arriving from 
each reactor will necessarily be quite low. This 
means both that the detector must be very large, 
appreciably larger than current proposals for RENO 
50 [8], and also that flux from multiple reactors 
must be used. Multiple reactors are available, es- 
pecially in places like Japan, Korea and China's 
Guangdong province where such experiments may 
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I K '.. 1: Simulated average values of RL + PV obtained 
from 100,000 neutrinos observed at a baseline of 58 km 
assuming the numerically interpolated reactor spectra 
from the 1980's [18-20]. The black dotted curve uses 
an unweighted Fourier transform and the red solid curve 
uses the weig ht exp(-0.04£; 2 /MeV 2 ). 



be built. However multiple reactors imply multi- 
ple baselines, and so vs will arrive at the detector 
with their 1-3 oscillations out of phase. Neutrinos 
from different reactors are not coherent with each 
other, wavefunctions are not added, but probabili- 
ties are added and this destroys the fine structure of 
the spectrum whose precise measurement is essential 
to a determination of the hierarchy. 

We now provide the most systematic analysis of 
this interference effect to date. Our goal is to illus- 
trate the effect of multiple baselines on the chance 
of success of the experiment and on the optimal lo- 
cation. We will consider two 18 GW thermal capac- 
ity idealized point reactor neutrino sources with flux 
given by the Gaussian form of Ref. [20], bearing in 
mind that a realistic flux with a weighted transform 
will lead to similar results. By defining an effective 
baseline difference, our results can be applied to con- 
figurations with many reactors. The reactor flux is 
normalized such that, including v oscillations, at 58 
km each 18 GW complex yields 25,000 I7 e in 3 years. 

In Fig. 2 we display the probability p of suc- 
cessfully determining the hierarchy with 3 and 6 
years of live time for various combinations of base- 
lines. The solid and dashed curves are analyses us- 
ing respectively RL + PV and a neural network op- 
timizing 24 coefficients corresponding to the 4 hi- 
erarchy indicators of Ref. [12] and to 6 weights: 

e -0.08B 2 e -0.02B 2 e -0.04(.E-3.6) 2 e ~E/S e -0.1(£-5.25) 2 

and e -(E-3-e) 3 /wo_ 

The first two panels are the 6 and 3 year probabil- 
ities of success for baseline differences between and 
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FIG. 2: Chance of success, p, after 6 years (first and third 
panels) and 3 years (second and fourth panels) of live 
time, with neutrinos from two 18 GW reactor complexes 
at distinct baselines. The horizontal axis is the shortest 
baseline, the color is the difference between baselines. 
On the first two panels these are (black), 500 m (red), 
1 km (blue), 1.5 km (purple) and 2 km (green). On 
the last two panels these are 5 km (black), 10 km (red), 
15 km (blue) and 20 km (purple). The solid curves use 
RL + PV and the dashed curves a neural network. 



4 



2 km. A 2 to 5 km baseline difference causes a lower 
chance of success. With these analysis methods, a 
95% chance of a successful hierarchy determination 
in 6 years is only possible if the baseline difference is 
appreciably less than 1 km, about the effective differ- 
ence for the BaiMianShi and Mudeungsan sites [21] 
for Daya Bay II and RENO 50. As RENO 50 uses 
neutrinos from only 18 GW of reactors, the times 
reported here must be doubled, in addition to the 
factor of four discussed earlier. 

The last two panels illustrate that, for baseline 
differences of 5 to 20 km, the far detector is a back- 
ground. As a result a larger baseline difference and 
a shorter baseline to the near reactor both increase 
the signal/background ratio and so p. These strong 
backgrounds reduce the optimal baseline. 

Interference and flux from more distant reactors 
is a big problem for differences above 500 meters, 
arising when a detector is not perpendicular to a long 
linear array like the Daya Bay - Ling Ao complex. 
Sites such as BaiYunZhang are perpendicular and 
so enjoy identical baselines, but use flux from only a 
single array. In particular Daya Bay II sites such as 
BaiMianShi are not competitive. Reno 50 will face 
similar problems, as the largest mountains are at a 
45 degree angle to the reactor array. 

The minimal baseline difference at potential Daya 
Bay II sites that use flux from multiple reactor com- 
plexes is about 500 meters, corresponding to the 
area around the Dongkeng site of Ref. [21] which 
uses flux from the Taishan and Yangjiang complexes. 
Thus we find that the best case probability of deter- 
mining the hierarchy is about 98%. This result de- 
pends upon the values of the neutrino mixing param- 
eters, in general we have found that a ler increase in 
sin 2 (#i 3 ) or [AMf^ or a ler decrease in |AMf 2 | can 
improve the hierarchy determination by O.lcr — 0.3cr. 

The disappearance channel at NO^A may provide 
a 1 — 2% determination of the atmospheric effective 
mass of Rcfs. [7, 22]. For 8 <~ n, this mass dif- 
fers from the high energy reactor effective mass (2) 
by about 1.5% and from the low energy reactor ef- 
fective mass lAMfil by nearly 3%. As the sign of 
these differences depends upon the hierarchy, NO^A 
disappearance data can improve the hierarchy deter- 
mination at Daya Bay II. However there is no such 
advantage if 5 ~ 0, as the atmospheric effective mass 
would be nearly equal to that of Eq. (2), which will 



measured more precisely at Daya Bay II. 
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